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Abstract 

The generahzed sine-Gordon (sG) equation was derived as an integrable generalization of 
the sG equation. In this paper, we develop a direct method for solving the generalized sG 
equation without recourse to the inverse scattering method. In particular, we construct 
multisoliton solutions in the form of parametric representation. We obtain a variety 
of solutions which include kinks, loop solitons and breathers. The properties of these 
solutions are investigated in detail. We find a novel type of solitons with a peculiar 
structure that the smaller soliton travels faster than the larger soliton. We also show that 
the short pulse equation describing the propagation of ultra-short pulses is reduced from 
the generalized sG equation in an appropriate scaling limit. Subsequently, the reduction 
to the sG equation is briefiy discussed. 

MSG: 35Q51, 37K10, 37K40 
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1. Introduction 

We consider the following generalized sine-Gordon (sG) equation 

Utx = (1 + ^dl) sin u, (1.1) 

where u — u{x,t) is a scalar- valued function, i/ is a real parameter, — d^/dx^ and 
the subscripts t and x appended to u denote partial differentiation. The generalized sG 
equation has been derived for the first time in [1] using bi-Hamiltonian methods. Quite 
recently, equation (1.1) with i/ < was shown to be a completely integrable partial 
differential equation (PDE) [2]. Indeed, constructing a Lax pair associated with it, the 
initial value problem of the equation was solved for decaying initial data. In the process, 
the Riemann-Hilbert formalism was developed to obtain eigenfunctions of the Lax pair. 
Soliton solutions are obtainable in principle, but their derivation needs a very complicated 
procedure. Although some qualitative features of traveling-wave solutions are discussed in 
a different context from the Riemann-Hilbert formalism, explicit expressions of solutions 
are not available as yet. 

The purpose of this paper is to obtain exact solutions of equation (1.1) with v < Q 
and discuss their properties. We consider real and nonperiodic solutions. In our analysis, 
we take v — —1 without loss of generahty. The exact method of solution used here is 
the bilinear transformation method which is a very powerful tool in obtaining special 
solutions of soliton equations [3, 4]. The method has wide applications ranging from 
continuous to discrete soliton equations. The central problem in the bilinear formalism 
is the construction of tau-functions which are introduced through dependent variable 
transformations . 

This paper is organized as follows. In section 2, we develop an exact method of 
solution. Specifically, we use a hodograph transformation to transform equation under 
consideration into a more tractable form. The transformed equation is further put into 
a system of bilinear equations by introducing appropriate dependent variable transfor- 
mations. We then construct explicit solutions of the bilinear equations by means of a 
standard procedure in the bilinear formalism. The multisoliton solutions are obtained in 
the form of parametric representation. In section 3, we describe properties of solutions. 
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First, we consider 1-soliton solutions which include kink and loop soliton solutions as well 
as a new type of multivalued functions. Throughout this paper we use the term "soli- 
ton" as a generic name of elementary solutions such as kink, loop soliton and breather 
solutions. A novel feature of regular kink solutions is found which has never been seen 
in the sG kinks. Next, the 2-soliton solutions are discussed. We address the kink- kink 
and kink-loop soliton solutions together with the 1-brcathcr solution. Last, we explore 
the general multisoliton solutions. Our particular concern is the multikink solution for 
which the large-time asymptotic is derived and the associated formula for the phase shift 
is obtained. A recipe for constructing the multibreather solutions is briefly described. As 
examples of the multisoliton solutions, we present a solution describing the interaction 
between a soliton and a breather as well as a 2-breather solution which are reduced from 
the 3 and 4-soliton solutions, respectively. In section 4, we point out a close relationship 
between the generalized sG equation and the short pulse equation which models the prop- 
agation of ultra-short optical pulses. We show that the generalized sG equation is reduced 
to the short pulse equation by taking an appropriate scaling limit. The parametric mul- 
tiloop soliton solution of the short pulse equation presented in [5] is reproduced from the 
corresponding one for the generalized sG equation. The similar limiting procedure is also 
applied to the formula for the phase shift. Subsequently, the reduction of the generalized 
sG equation to the sG equation is discussed shortly. Section 5 is devoted to conclusion. 
In the appendix A, we show that the tau-functions for the multisoliton solutions obtained 
in section 3 satisfy a system of bilinear equations. The proof is carried out by means 
of an elementary method using various formulas for determinants. In the appendix B, 
we derive the 1-soliton solutions by an alternative method and demonstrate that they 
reproduce the corresponding 1-soliton solutions obtained in section 3. 

2. Exact method of solution 

2.1. Hodograph transformation 

We introduce the new dependent variable r in accordance with the relation 

r^^l + ul, (r>0), (2.1) 
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to transform equation (1.1) with v — —1 into the form 

rt + (r cos u)^ = 0. (2.2) 
We then define the hodograph transformation {x, t) {y, r) by 

dy — rdx — r cos u dt, dr — dt. (2.3) 

The X and t derivatives are then rewritten in terms of the y and r derivatives as 

d d d d d 

dx dy^ dt dr dy' 

With the new variables y and r, (2.1) and (2.2) are recast into the form 

= 1 + r\J, (2.5) 

(^i^ -{cosu)y = 0, (2.6) 
respectively. Further reduction is possible if one defines the variable by 

Uy = sin (p, (p = (p{y,T), (^-- < < -, mod 27rj . (2.7) 

It follows from (2.5) and (2.7) that 

- = cos 0. (2.8) 
r 

Substituting (2.7) and (2.8) into equation (2.6), we find 

0^ = sin u. (2.9) 

If we eliminate the variable (p from (2.7) and (2.9), we obtain a single PDE for u 

^^si„«. (2.10) 

Similarly, ehmination of the variable u gives a single PDE for (f) 



^1^ 



sin 0. (2.11) 



By inverting the relationship (2.4) and using (2.8), equation that determines the inverse 
mapping (y, r) — > (x, t) is found to be governed by the system of hnear PDEs for x — 

x{y,T) 



Note that the integrabihty of the system of equations (2.12) is assured by (2.7) and (2.9). 

A sequence of transformations described above are almost the same as those employed 
for solving the short pulse equation [5]. The underlying idea is to transform the orignal 
equation to the (possibly) integrable equation. In the case of the short pulse equation, the 
transformrd equation is the sG equation whereas in the present case, the corresponding 
equations are (2.10) and (2.11). The soliton solutions of the latter equations will be 
construcred here for the first time. Given u and 0, the most difficult problem is how to 
integrate the system of equations (2.12). This becomes the core part of the present paper 
and will be resolved by Theorem 2.1. 

2.2. Bilinear formalism 

Here, we develop a method for solving a system of PDEs (2.7) and (2.9). We use the 
bilinear transformation method [3, 4]. Let a and a' be solutions of the sG equation 



Xy — cos (f), 



(2.12a) 



x-j- — cos u. 



(2.126) 



— sm (7, 



a = cr{y,T), 



(2.13a) 



a = sm c7 , 



(7 — 



(2.136) 




(2.14a) 



(2.146) 



In terms of a and a', equations (2.7) and (2.9) can be written as 




(2.15a) 



(2.156) 



5 



It should be remarked that the system of PDEs (2.15) constitutes a Backlund transfor- 
mation of the sG equation with a Backlund parameter taken to be 1 [6]. The real- valued 
solutions of the sG equations (2.13) can be put into the form [7-10] 

(7 = 2ilny, (2.16a) 

(7' = 2iln^, (2.166) 

9 

where /* and g* denote the complex conjugate of / and g, respectively. The tau-functions 
/ and g play a central role in our analysis. They are fundamental quantities in constructing 
solutions. For sohton solutions, they satisfy the following bihnear equations [7]: 

D,Dyf.f=^{f-r'), (2.17a) 

D,Dyg-g^^{g'-g*'), (2.176) 
where the bilinear operators and Dy are defined by 

LT^.V -g-idr- dr^r (dy - dy,f f(r, y)g(r', y') \r'=r,y'=y, (m, n = 0, 1, 2, ...). 

(2.18) 

Now, we seek solutions of equations (2.7) and (2.9) of the form 

p* 

u^ihi—, (2.19a) 
F 

= iln— , (2.196) 

where F and G are new tau-functions. It turns out from (2.14), (2.16) and (2.19) that 

2ii = (j + fj' = 2i In^- = 2i In— , (2.20a) 

20 = cT-a' = 2iln4^ = 2iln^. (2.206) 
fg* G 

By taking into account (2.20), we may assume the following relations among the tau- 
functions f,g,F and G: 

F = fg, (2.21a) 
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G^fg*. 

The above expressions lead to an important relation 



(2.21&) 



F*F = G*G. (2.22) 

Substituting (2.19) into equations (2.7) and (2.9) and using (2.22), we see that F and G 
satisfy a system of bilinear equations 

DyF* ■ F ^ ~{G^ -G*^), (2.23a) 

DrG* - G = --{F^ - F*^). (2.236) 

Thus, the problem under consideration is reduced to obtain solutions of equations (2.23) 
subjected to the condition (2.22). After some trials, however, we found that this proce- 
dure for constructing solutions is difficult to perform. Hence, we employ an alternative 
approach. To begin with, we impose the following bilinear equations for / and g which 
turn out to be the starting point in our analysis 

Dyf-g*^l{fg*-rg), (2.24a) 
Drf-g^lifg-fg*). (2.246) 

With (2.24) at hand, the following proposition holds: 

Proposition 2.1. /// and g satisfy the bilinear equations (2.24), then the tau-functions 
F and G defined by (2.21) satisfy the bilinear equations (2.23). 

Proof. First, we prove (2.23a). We substitute (2.21a) into the left-hand side of (2.23a) 
and rewrite it in terms of bilinear operator to obtain 

DyF* ■ F = -{DJ ■ g*)f*g + {DJ* ■ g)fg\ (2.25) 

By virtue of (2.24a), the right-hand side of (2.25) becomes —{l/2){{fg*Y — {f*gY} which 
is equal to the right-hand side of (2.23a) by (2.21b). The proof of (2.23b) can be done in 
the same way by using (2.24b). □ 

2.3. Parametric representation 
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We demonstrate that the solution of equation (1.1) with u — —1 admits a parametric 
representation. The following relation is crucial to integrate (2.12): 
Proposition 2.2. cos is expressed in terms of f and g as 

cos0 = l+ (^In^^ . (2.26) 
Proof. Using (2.24a), one obtains 

.ll] Dyf-9* Dyr-9 

f*f)y fg* f*9 

2 f*fg*g ^ ' 

On the other hand, it follows from (2.19b) and (2.21b) that 

, I f G G 
cos = - — + 



2\G* G 

1 {fg*f + {rgf 



(2.28) 



2 f*fg*g 

The relation (2.26) follows immediately by comparing (2.27) and (2.28). □ 
Integrating (2.12a) coupled with (2.26) by we obtain the expression of x 

x^y + ln^ + dir), (2.29) 

where d is an integration constant which depends generally on r. The expression (2.29) 
now leads to our main result: 

Theorem 2.1. The real-valued solution of equation (1-1) with v — —1 can he written by 
the parametric representation 

uiy,T)^iln^, (2.30a) 

x{y,T)^y + T + lnj^+yo, (2.306) 

where the tau-functions f and g satisfy both (2.17) and (2.24) simultaneously and yo is 
an arbitrary constant independent of y and r. 

Proof. The expression (2.30a) for -u is a consequence of (2.19a) and (2.21a). To prove 
(2.30b), we substitute (2.29) into (2.12b) and obtain the relation 

cosu^(\n^^ +d!{T). (2.31) 
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The left-hand side of (2.31) can be expressed by / and g in view of (2.30a) whereas the 
right-hand side is modified by using (2.24b). After a few calculations, we find that most 
terms are cancelled, leaving the equation ^'(t) = 1. Integrating this equation, one obtains 
d{T) = T + Ho, which, substituted into (2.29), gives the expression (2.30b) for x. □ 
The parametric solution (2.30) would produce in general a multi-valued function as 
happned in the case of the short pulse equation [5] . To derive a criterion for single- valued 
functions, we calculate with use of (2.7) and (2.8) and obtain 



— rUy — tan (p. 



(2.32) 



Thus, if the inequality — 7r/2 < < 7r/2 (mod tt) holds for all y and r, then u becomes a 
regular function of x and t. By virtue of the identity 



as well as the relation (2.21b), the above condition for regularity can be written as 



(2.33) 



Mfg*) 



(2.34) 



In the case of 1-soliton solutions discussed in the next section, the above condition is 
found explicitly in terms of the parameters characterizing solutions. 

2.4- Multisoliton solutions 

The last step in constructing solutions is to find the tau-functions / and g for the sG 
equation which satisfy simultaneously the bilinear equations (2.24). The following theorem 
establishes this purpose: 

Theorem 2.2. The tau-functions f and g given below satisfy both the bilinear forms 
(2.17) of the gG equation and the bilinear equations (2.24) 

N 



/ = ^ exp 



.3=1 l<j<k<N 



5 = exp 

("=0,1 



,j=l l<j<k<N 



(2.35a) 



(2.356) 



where 



= PjV + -r + Cjo, U = 1, 2, N), 



Pj 



(2.36a) 




Pj - Pk 
Pj + Pk 



)■■ 

-Pj 
+ Pj 



{j,k^l,2,...,N;j^k), 



(i = 1,2, ...,7V). 



(2.366) 



(2.36c) 



Here, Pj and C,jo are arbitrary real parameters satisfying the conditions pj ^ pk for j ^ k 
and N is an arbitrary positive integer. The notation X]^=o i ^"^P^^es the summation over 
all possible combination of jii = ^2 = 0, 1, jijq = 0, 1. 

Proof. It has been shown that / and g given by (2.35a) and (2.35b) satisfy the bihnear 
equations (2.17a) and (2.17b), respectively [7]. Thus, it is sufficient to prove that they 
satisfy the bihnear equations (2.24). The proof is carried out by a lengthy calculation 
using various formulas for determinants. It will be summarized in the appendix A. □ 

2.5. Remark 

The tau-functions / and g given by (2.35) yield real-valued solutions since all the param- 
eters Pj and {j — 1, 2, N) are chosen to be real numbers. If one looks for breather 
solutions, for example, one needs to introduce complex parameters (see sections 3.2.3, 
3.3.2 and 3.3.3). Even in this case, however, the analysis developed here can be applied 
as well without making essential changes. Actually, we may use the tau-functions /' and 
g' instead of /* and g*, respectively where /' and g' are obtained simply from / and 
g by replacing i by — i, but all the parameters in the tau-functions are assumed to be 
complex numbers. The solutions of the sG equations (2.13a) and (2.13b) can be written 
as a — 2i ln(/'//) and a' — 2i \n{g'/g), respectively, where the tau-functions f,f,g and 
g' satisfy the following systems of bilinear equations 



DrDyf ■ f = - f"), D^Dyf . f = ^if" 
D^Dyg .g = hg^-g''), D^Dyg' ■g' = l {g" 




n 



(2.37a) 



(2.376) 



The bilinear equations corresponding to (2.24) are then given by 



Dyf-g' = l{fg'-f'g), 



Drf-9 = ^{f9-fg'), 



(2.38a)) 
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Dyf ■ 9 = \{fg - fg'), DJ' ■ g' = i(/V - fg). (2.386)) 

The expressions corresponding to (2. 19)- (2.23) are obtained if one replaces the asterisk 
appended to the tau-functions by the prime. Under these modifications, the real-valued 
solutions are produced if one imposes the conditions /' = /* and g' — g*. See also an 
analogous problem associated with breather solutions of the short pulse equation [5] . 

3. Properties of solutions 

The parametric representation (2.30) of the solution with the tau functions given by (2.35) 
exhibit a variety of soliton solutions of equation (1.1) with u = —1. As exemplified here, 
solutions include both single-valued and multi-valued kinks, loop solitons and breathers. 

3.1. 1-soliton solutions 

The tau-functions for the 1-soliton solutions are given by (2.35) with N = 1: 

/ = l + ie«S = Pil/ + - + Cio, (3.1a) 

Pi 

g = l + is,e^\ s, = l^. (3.16) 

1 +Pi 

The parameters pi and are related to the amplitude and phase of the soliton, re- 
spectively and ^1 is a phase variable characterizing the parametric representation of the 
solution. It follows from (2.21) and (3.1) that 

F = 1 - sie^^i + i (1 + si)e^i , (3.2a) 

G = 1 + sie^^i + i (1 - si)e«i . (3.26) 
Using (2.30) and (3.1), the parametric representation of the solution is written in the form 

u — 2 tan~^(sinh — pi cosh ^i) -|- tt, (3.3a) 

X = y + T + In f - tanh + yo, (3.36) 

where we have imposed the boundary condition «(— oo,t) = 0. Note that if u solves 
equation (1.1), then so do the functions ±u + 27in (n : integer). To describe solutions of 
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traveling-wave type like 1-soliton solutions, it is convenient to parameterize the solutions 
in terms of the single variable ^i. For this purpose, we rewrite (3.3b) as 



X = x + Cit + Xo^ h In — - — tanh + yo, 

pi V(l+Pi) (1+Pi) / 



(3.4a) 



where ci is the velocity of the soliton given by 



ci = ^ - 1, (3.46) 
Pi 

and xq = iw/pi- Observing the motion of the soliton in the original (x, t) coordinate 
system, it travels to the left at the constant velocity Ci for pi < 1 and to the right for 
Pi > 1. In the critical case pi = 1 for which ci = 0, the soliton remains stationary. The 
profile of the soliton changes drastically depending on values of the parameter pi. The 
singular nature of the solution can be extracted conveniently from the information on the 
gradient of u with respect to X. A calculation using (3.3a) and (3.4) gives 

1 I- Pi tanh ^1 

ux = 7—2 : — 2- (3-5) 

pi cosh ^1 tanh^ ^1 - ^ tanh 6 + ^ 

Let us now analyze the solution (3.3). Using (3.1), the condition (2.34) for single- 
valued function becomes 

-1 < ■ , , < 1. (3.6) 

cosh ^1 — pi smh 

In the following analysis, we assume pi > without loss of generality. If < pi < 1, then 
(3.6) leads to the inequality 



1 l-v^ 

tanh^ ^1 tanh H > 0. (3.7) 

Pi 2p^ 



One can see that the inequality (3.7) always holds for < pi < In this case, ux from 
(3.5) becomes finite for arbitrary values of On the other hand, if ^ < pi < 1, then 
the solution exhibits singularities at two points X = X± = X{^±) where 



^± = tanh 



(3.^ 
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If Pi > 1, the inequality (3.6) breaks down, giving rise to two singular points whose 
positions are the same as X±. Unlike the second case, however, ux becomes zero at 
X — Xo = X{^o) with ^0 = tanh~^(l/]9i), as readily noticed from (3.5). The solution for 
the case pi = 1 exhibits a peculiar behavior, which deserves a separate study. In view of 

these observations, the solutions can be classified to four types according to values of pi, 
or equivalently ci by (3.4b), which we shall now investigate in detail. 

Type 1. Regular kink: Ci > 1 (0 < pi < 

Figure 1 shows a typical profile of as a function ol X. It exhibits a profile of a 27r-kink 
similar to the kink solution of the sG equation [6]. The propagation characteristic is, 
however, different from that of the sG kink. To clarify this point, we rewrite ux from 
(3.5) as 

ux = cosh(e.-dO 

yr^cosh^(6-rfo-4. 

where di = (1/2) ln[(l +pi)/(l — pi)]. It turns out from (3.9) that the amplitude A of ux 
is given by 

A^'J^. (3.10) 

Eliminating pi from (3.4b) and (3.10), we find the dependence of the amplitude on the 
velocity 

ci = ^(a^ + 2 + 2VA^ + l) . (3.11) 

The relation (3.11) indicates that ci is a monotonically decreasing function oi A. In other 
words, the smaller soliton travels faster than the larger sohton. This pecuhar feature of 
the solution has never been observed in the behavior of the sG kink solutions. The broken 
line in figure 1 plots the profile of v = ux obtained from the kink solution depicted in the 
same figure. It represents a soliton. 

Figure 1 

Type 2. Singular kink: < ci < 1 (^ < pi < 1) 

The profile of the solution of type 2 is a 27r kink, but in the interval X_ < X < X_|_, it 
becomes a three- valued function. Figure 2 shows a typical profile of m as a function of X. 
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Figure 2 

Type 3. Loop soliton: — 1 < Ci < {pi > 1) 

When compared with solutions of types 1 and 2, the solution for pi > 1 exhibits a 
different behavior. Indeed, we see from (3.3), and (3.5) that m(±oo) = 0,Ux{Xq) = 0, 
respectively. In addition, the solution has two singular points at X — X± such that that 
Xq — (X+ + XJ)/2. Figure 3 shows a typical profile of as a function ol X. It represents 
a loop soliton and has a symmetrical profile with respect to a straight fine X — Xq. 

Figure 3 

Type 4- Stationary solution: Ci — 1 {pi — 0) 

For this special value of the parameter pi, the parametric solution (3.3a) and (3.4a) takes 
the form 

tan-^ e~^i + tt, (3.12a) 

X = - \n{e^' +e-^'). (3.126) 

Remarkably, one can eliminate the variable from this expression to give an explicit 
solution 

sinM = 2e^, (-oo < X < - In 2). (3.13) 

The profile of u is illustrated in figure 4 as a function of X. Since ci = 0, one has 
X — X + Xq by (3.4a) and the solution becomes time-independent. Keeping this fact in 
mind, we can derive the solution (3.13) directly from the stationary version of equation 
(1.1) with i/ = -1 



l-^)-» = 0. (3.14) 

Integrating equation (3.14) under the boundary condition u{—oo) — (modvr), we recover 
(3.13). 

Figure 4 

The 1-soliton solutions presented above are of fundamental importance in constructing 
general multisoliton solutions. In fact, the latter solutions will be shown to consist of any 
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combination of the former solutions. In appendix B, we shall derive the 1-soliton solutions 
by means of an elementary method. 

3.2. 2-soliton solutions 

The tau-functions for the 2-soliton solutions read from (2.35) with N — 2 in the form 



where 5 = (pi -p2)V(Pi ^Pif and = (1 -pi)/(l +pj) {i = 1,2). 

The parametric solution (2.30) represents various types of solutions describing the in- 
teraction of two solitons. Here, soliton means one of four types of solutions given in section 
3.1. In addition, the solution exhibits a breather solution when the two parameters pi 

and p2 appear as a complex conjugate pair. Although a number of solutions yield accord- 
ing to the combination of elementary solutions, we present three types of solutions, i.e., 
kink-kink (or two-kink) solution and kink-loop soliton solution and a breather solution. 

3.2.1. Kink-kink solution 

The solid line in figure 5 a-c exhibits the profile of solution u for three different times. 
The solution represents the 47r-kink. In the same figure, we depict v = Uxhj the broken 
hue, showing that it represents the interaction of two sohtons. The characteristic of the 
interaction process of two solitons is seen to be quite different from that of the sG two 
solitons. Indeed, as evidenced from figure 5, a smaller soliton overtakes, interacts and 
emerges ahead of a larger sohton. This reflects the fact that the velocity of each sohton is 
a monotonically decreasing function of its amplitude (see (3.11)). After the interaction, 
both solitons suffer phase shifts. The general formula for the phase shift arising from the 
interaction of N solitons will be given by (3.26) below. In particular, for = 2, it reads 



/ = l + i(e^i +e^^) - (5e^l+^^ 



(3.15a) 



g = 1 + i(sie^^ + 526^^) - 5siS2e' 



(3.156) 




(3.16a) 




(3.166) 
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In the present example, formula (3.16) yields Ai = 4.55 and A2 = —2.42. A careful 
inspection of (3.16) reveals that Ai > andA2 < for arbitrary values of pi and p2 



has moved forward and the large soliton backward relative to the positions they would 
have reached if both solitons had moved at constant velocities throughout the interaction 
process. The novel feature of the 2-soliton solution described above will appear here for 
the first time. 

Figure 5 a-c 

3.2.2. Kink-loop soliton solution 

The next example is a solution u representing the interaction of a 27r-kink and a loop 
sohton. See figure 6 a-c. Since the kink propagates to the left and the loop soliton to the 
right, the solution describes the head-on collision unlike the first example which exhibits 
an overtaking collision. 



3.2.3. Breather solution 

The breather solution can be interpreted as a bound state composed of a kink and antikink 
pair in the sG model [6]. It has a localized structure which oscillates with time and decays 
exponentially in space. In the generalized sG equation, the similar breather solutions to 
the sG breathers will be shown to exist. The procedure for constructing breather solutions 
is the same as that has been used for the short pulse equation [5]. To be more specific. 



satisfying the inequality < pi < p2 < l/v^- This implies that the small soliton 



Figure 6 a-c 



let 



Pi = a + ih, 



P2 = a-ib = pI, 



(a > 0, b> 0), 



(3.17a) 



^10 = A i//. 



60 = A-i/x = ^i, 



(A, II : real). 



(3.176) 



Then, / and g from (3.5) become 




(3.18a) 




(3.186) 
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where 



6 = ^ + iX, 



(3.18c) 




(3.18(i) 



(S.lSe) 




(3.18/) 




(3.19a) 



(3.196) 



The regular solution is obtainable if the inequality (2.34) holds with / and g being given 
by (3.19). An inspection reveals that if a/b is sufficiently small compared to 1, then the 
solution would exhibit no singularities. However, it is not easy to extract the condition 
for the regularity from (2.34) when compared with the corresponding problem for the 
breather solution of the short pulse equation [5]. We leave it to a future work. Instead, 
we present a regular solution by a numerical example. In figure 7a-c, the profile of v = 
is depicted for three different times. We can observe that the breather propagates to the 
left while changing its profile. The propagation characteristic of the breather is similar to 
that of the short pulse equation [5] . 



3.3. N-soliton solutions 

The solutions including an arbitrary number of solitons can be constracted from (2.30) 
with the tau-functions (2.35). There exist a variety of solutions which are composed of 
any combination of 1-soliton solutions presented in section 3.1. Here, we address the N- 
kink solutions and M{— 2N) breather solutions. For the former solutions, we investigate 
the asymptotic behavior of solutions for large time and derive the formulas for the phase 



Figure 7a-c 
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shift while for the latter ones, we provide a recipe for constructing M breather solution 
from the A'"-soliton solution. As examples, we present a solution describing the interaction 
between a soliton and a breather as well as a two-breather solution. 



3.3.1. N-kink solution 

Let the velocity of the jth kink be Cj = — 1, (0 < pj < and order the 

magnitude of the velocity of each kink as ci > C2 > ... > ctv- We observe the interaction 
of kinks in a moving frame with a constant velocity c„. We take the limit t — > — oo 
with the phase variable ^„ being fixed. We then find that / and g have the following 
leading-order asymptotics 

TV 1 

f^S^exp J2 (^i + iO (l + ie^"^'"), (3.20a) 



U=n+1 



5f ~ 5n exp 



N 1 

E (0-2ci, + |i) (l-ie^"-'^"-^^^-'). 



J=n+1 



where 



n 



Pj - Pk 
Pj + Pk 



(3.206) 

(3.20c) 
(3.20(i) 



n+l<j<k<N 

If we substitute (3.20) into (2.30), we obtain the asymptotic form of u and x: 



u 2 tan ^ 



sinh(^„ - dn + Si ^) 



cosh dr, 



N 



l+p„tanh(^„-rf„ + 5i ^) , , , 

X ~ y T - In — — ; ^ - 4 2^ dj- 2d„ + yg. 



(3.21a) 



(3.216) 



1-Pn tanh(^„ -dn + 6k j j^^i+i 
Note that we have used equivalent but different expressions for u and x from those given 
by (3.3). 



As i — > -|-oo, the expressions corresponding to (3.21) are given by 

sinh(^„ - dn + Sii'^] 



w ~ 2 tan ^ 



cosh dr, 



+ 7r. 



(3.22a) 
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1 +y;„,taiili(4„ - r/„ + . , ^, , 

x^y + T -In -4 >rfj -24 + 7/0- (3.226) 

1-Pn tanh(^„ -dn + Sk ') j=i 

with 

e'^|:^(^^)^ (3.20 

Let be the center position of the nth kink in the (x, t) coordinate system. It simply 

stems from the relations ~ dn + Slt^ = by invoking (3.21a) and (3.22a). Thus, as 
t — oo 

1 _ ^ 

Xc + cj + x„o —{dn- Sl~^) - 4 ^ dj- 2d„ + yo, (3.23) 

^" j=n+l 

where Xno — Cno/Pn- As t ^ +oo, on the other hand, the corresponding expression turns 
out to be 

^ n—l 

Xc + Cnt + Xno (^n - 5^+^) - 4 ^ d,- - 2dn + yo, (3.24) 

If we take into account the fact that all kinks propagate to the left, we can define the 
phase shift of the nth kink as 

A„ = Xc{t -oo) - Xc{t +oo). (3.25) 

Using (2.36c), (3.20c), (3.22c), (3.23) and (3.24), we find that 

A. = -\yin - y In ^-^) 

Pn [jr'i \Pn+PjJ \Pn+PjJ 

n-1 / -, . \ 2 A'' 



(3.26) 

j=l V- -f-// j=n+l ^ ' "" ' 

The first term on the right-hand side of (3.26) coincides with the formula for the phase 
shift arising from the interaction of N kinks of the sG equation [7, 8, 10] whereas the 
second and third terms appear as a consequence of the coordinate transformation (2.3). 

3.3.2. M-breather solution 

The construction of the M-breather solution can be done following the similar procedure 
to that for the 1-breather solution developed in section 3.2.3. Here, M(= 2N) is a 
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positive even number. To proceed, we specify the parameters in (2.35) and (2.36) for the 
tau-functions / and g as 

P2j-i = p*2j = aj + ibj, aj > 0, bj > 0, (j = 1, 2, M), (3.27a) 

6,-1,0 = C2j,o = A,- + il^j, {j = 1, 2, M), (3.276) 

where Uj and bj are positive parameters and Xj and //j are real parameters, respectively. 
Then, the phase variables ^j-i and are written as 

6,-1 = 0j + iXj, U = 1, 2, M), (3.28a) 

C2j^ej-ixj, (j = l,2,...,M), (3.286) 

with 

= a^iy + CjT) + A„ (j = 1, 2, M), (3.28c) 

X,- = 6,(y - c,t) + /x^, (j = 1, 2, M), (3.28ci) 

^j-^rn^^ (i = l'2,...,M). (3.28e) 

^ 

The parametric solution (2.30) with (3.27) and (3.28) describes multiple collisions of M 
regular breathers provided that certain condition is imposed on the parameters aj and 
bj{j — 1,2,...,M). Although it will be a difficult task to derive the condition for the 
regularity through the inequality (2.34), numerical examples confirm the existence of 
regular multibreather solutions. See also an example of the 2-breather solution of the 
short pulse equation [5]. The asymptotic analysis for the M-breather solution can be 
performed following the procedure for the A^-kink solution, showing that the M-breather 
solution splits into M single breathers as t — )■ ±oo. The resulting asymptotic form is, 
however, too complicated to write down and hence we omit the detail. One can refer to 
the similar analysis to that for the M-breather solution of the short pulse equation [5] . 

3.3.3 Soliton-breather solution 

We take a 3-soliton solution with parameters pj and 6j (i = 1? 2, 3). If one impose the 
conditions p2 — Pi, 62 — Coi ^ already specified for the breather solution (see section 
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3.2.3) and < ps < ^03(:real) for the regular kink solution, then the expression 

oi V = Ux would represent a solution describing the interaction between a soliton and a 
breather. We choose pi,p2,Cio and ^20 as those given by (3.17). Then, the tau-functions 
/ and g from (2.35) become 

/ = l+i(e^i+e^i*+e^3)+Q e^^+^i-^ige^i+^^^^^gC^+^^+i (^^^ 5l35^3e^^+^*+^^ (3.29a) 
g = l + i(sie«i + sle^i + 536^^) + sisle^'+^'i - Si^siSse^'^^^ - 5^351^36^1+^^ 

^'{l) ^l35^3Sls^s3e«^+«+«^ (3.296) 

where 

l-a-i6 ^ I-P3 /a-p3 + i6y ro on \ 

si = 1 , , = S3 = 7— — , 5i3 = — —77 = ()23. 3.29c) 

l + a + 16 I+P3 \a + p3 + i6/ 

Figure 8a-c shows a profile of f = m^, for three different times. Wc see that as time 
evolves the soliton overtakes the breather whereby it suffers a phase shift. An asymptotic 
analysis using the tau-functions (3.29) leads to the formula for the phase shift of the 
soliton, which we denote A. Actually, one has for p| < + 6^ 



2 (P3 + ar + , , (1 + ay + ,^ , 

^ = — ^^7 ^^~w~2 1n- ^—jo- 3.30a) 

P3 (P3 - a) +0 (1 - a) + ^ 



and for a^ + 6^ < pi 



2 , (p3 + a)2 + 62 (i + a)2 + 52 

^ = 1^ 7^ + 2 1^ 71 kn^- 3.306 

P3 (P3 - «) +0"' (1 - ay + 6^ 

In the present example, formula (3.30a) gives A = 3.07. 

Figure 8a-c 

3.3.4 Breather-breather solution 

The breather-breather (or 2-breather) solution is reduced from a 4-soliton solution fol- 
lowing the procedure described in section 3.3.2. Figure 9a-c shows a profile oi v = for 
three different times. It represents a typical feature common to the interaction of solitons, 
i.e., each breather recovers its profile after collision. 
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Figure 9a-c 



4. Reduction to the short pulse and sG equations 

The short pulse equation was proposed as a model nonlinear equation describing the 
propagation of ultra-short optical pulses in nonlinear media [11]. It may be written in an 
appropriate dimensionless form as 

utx^u + hu^)^^, (4.1) 



where u — u{x, t) represents the magnitude of the electric field. Here, we demonstrate that 
the generalized sG equation is reduced to the short pulse equation by taking an appropriate 
scaling limit combined with a coordinate transformation. The A'"-soliton solution of the 
short pulse equation as well as the formula of the phase shift can be derived from those of 
the generalized sG equation. The reduction to the SG equation is shown to be established 
as well by another scaling limit. 

4.I. Reduction to the short pulse equation 

Let us first introduce new variable it, x and i according to the relations 

- u _ 1 , , - , ^ . 

u^-, x^-[x-t), t^et, (4.2) 

where e is a small parameter and the quantities with bar are assumed to be order 1. 
Rewriting the derivatives in terms of the new variables t and x and expanding sin eu in 
an infinite series with respect to e, we can develop equation (1.1) with u — —1 to 



e ( - \uxx-\ ^ eu - —u^ + ... - ( eu - —u^ + —u^ + ...]. (4.3a) 
\ tx ^2 6 e2 V 6 120 ' ^ ^ 



Note that the terms of order e ^ are canceled. Thus, we have 

eutx = e(^u+ ^{u%^)^ + 0{e^). (4.36) 

If we divide both sides of (4.3b) by e and then take the limit e ^ 0, we arrive at the short 
pulse equation (4.1) written by the new variables. It is noteworthy that the key relation 
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(2.1) has been used to transform the short pulse equation into the sG equation and it is 
invariant under the scahng (4.2). The similar scahng variables to (4.2) have been used 
to derive the short-wave models of the Camassa-Holm and Degasperis-Procesi equations 
[12]. 

4.1.1. Scaling limit of the N-soliton solution 

In order to perform the scaling limit of the A^-soliton solution, we find it appropriate to 
employ the following new variables in addition to the variables defined by (4.2) 



- y - yo 
y = - yo = — , r = er, 
e e 



Pj = ^Pj: Oo = Oo, (j = l,2,...,iV). (4.4) 
If one rewrites the tau- function / from (2.35a) in terms of these variables, one simply has 



f = f = ^ exp 

M=0,1 



N 



with 



r 



l<j<k<N 



= PjV + - + ^jo, (j = l,2,...,iV), 
Pj 



Pj - Pk 



{j,k = l,2,...,N;jj^k). 



(4.5a) 



(4.56) 



(4.5c) 



^Pj + Pk , 

In proceeding to the limiting procedure for the tau-function g, we need to retain terms 
up to order e. Thus, the expansion 



Af \ / N 



(4.6) 

as well as the scaled variables (4.2) and (4.4) are substituted into (2.35b) to derive the 
expansion of g 



N 



^ ^ ^ j 1 _ -f- I exp 



■ N 

.j=l l<j<k<N 



fJ'jfJ'klfjk 



= f-2e/? + 0(e^). 



(4.7) 
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Insertion of (4.2), (4.5) and (4.7) into (2.30a) gives 



In the limit of e ^ 0, (4.8) leads to the scaling limit of u 



M = 2i ( In 4 ) ■ (4-9) 



Applying the similar procedure to (2.30b), we find 

ex^ey- 2e(ln /*/), + eyo + 0{e^). (4.10) 

It follows from (4.10) that 

x^y-2{\nrf)r + yo. (4.11) 

The expressions (4.9) and (4.11) coincide with the parametric representation of the A^- 
soliton solution of the short pulse equation [5]. 

4.1.2. Scaling limit of the phase shift 

The scahng limit of the formula (3.26) for the phase shift can be derived easily. Indeed, if 
we define the new variable for the phase shift by A„ = A„/ e and substitute this expression 
and the scaled variable pj = epj from (4.4) into (3.26), we find, after taking the limit e — > 0, 
that 

Pn [j^^ \Pn+PjJ \Pn+Pj. 
n-l . ^ A 

+ E (n^l,2,...,N). (4.12) 

j=l j=n+l 

This expression is just the corresponding formula for the short pulse equation [5]. Note 
that if all pj are real parameters such that pj 7^ pk for j k, then (4.12) gives the formula 
for the phase shift resulting from the overtaking coUisions of N loop solitons. 

4.2. Reduction to the sG equation 
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The reduction to the sG equation is rather straightforward compared to the previous one 
for the short pulse equation. It turns out that the appropriate scaled variables are given 

by 

- t _ T 

u — u, X = ex, y = ey, t — -, t — 

e e 

Pj = ^, ^,-o = U (i = 1,2, ...,7V). (4.13) 

In terms of the variables (4.13), we can recast equation (1.1) to the sG equation %x- = sin-u 
in the limit of e — ?■ 0. The parametric solution (2.30) reduces to the usual form of the 
iV-soliton solution of the sG equation i.e., u{x,t) = 2i ln(/*//) where / is given by (4.5) 
with the identification y = x,f = t. The phase shift is scaled by A„ = eA„. It is given 
by the first term on the right-hand side of (4.12), reproducing the well-known formula 
derived by the asymptotic analysis of the A^-sohton solution of the sG equation [7, 8, 10]. 

5. Conclusion 

A direct approach employed in this paper constructs various types of soliton solutions 
such as single- and multi-valued kinks, loop solitons and breathers. These elementary 
solutions are combined to produce a variety of multisoliton solutions. As examples, we 
presented a solution describing the interaction between a soliton and a breather as well 
as a 2-breather solution. As far as solutions are concerned, one can observe that the 
generalized sG equation has a rich structure compared with that of the sG equation. We 
also demonstrated that the generalized sG equation is reduced to both the short pulse and 
sG equations in appropriate scaling limits. Another interesting issue is the generalized sG 
equation (1.1) with v — 1 for which the method of solution still remains open. One may 
apply a sequence of nonlinear transformation similar to that used here to obtain solutions. 
Another direction to be worth investigating is the periodic problem. The exact method of 
solution used here will work well for constructing periodic solutions of equation (1.1). See 
[13, 14] for periodic solutions of the short pulse equatoin. These problem will be pursued 
in future works. 
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Appendix A. Proof of Theorem 2.2 



In this appendix, we show that the tau-functions / and g given respectively by (2.35a) 
and (2.35b) satisfy the bihnear equations (2.24). First, we rewite / and g in terms of 
determinants. For the purpose, we use the formula [15] 



exp 

iu=0,l 



N 



.3=1 l<j<k<N 



Ajv det ie'^'5jk + 



2P, 



Pj + Pk/ l<j,k<N 



where 



N 



= ^j+Yl ^J*^' Aat = exp I - jjk ] , 

l<j<k<N 



(A.l) 



{A.2) 



k=l 



and 6jk is Kronecker's delta. Since numerical factors multiplied by / and g have no effects 
on the proof of (2.24), wc use the determinantal expression given by the right-hand of 
(A.l) instead of the finite sum. Furthermore, we shift the phase factor ^jq by — ^^=i 7jjk 
SO that Q — ^j. Consequently, we can express / and g by the following determinants: 



/ = det A = 1^1, A= {ajk)i<j,k<N, ajk = ie^'^jfc + 

X — T) ' 

g ^ det B = \B\, B = {bjk)i<j,k<N, bjk = i ^ e^'Sjk + 

J- ~r 



2pj 



Pj + Pk ' 
2pj 



{A.3) 
{AA) 



+ Pj Pj + Pk 

For later convenience, we introduce some notations as well as formulas for determi- 
nants. Matrices and cofactors associated with any N x N matrix A — {ajk)i<j,k<N are 
defined as follows : 



A{a;h) 



A{a., b; c, d) — 



^ ail ■ 


■ diN 


bl\ 




ajsfi . 


■ OjvAT 


bN 




\ ai . 


■ aN 






(ail ■ 


■ CllN 


Cl 


di\ 


dm ■ 


■ dNN 


cn 




ai . 


Oat 








\bi . 


• bN 








(A.5) 



^jk 



d\A[^ 
dajk ' 



(A.6) 
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Here, Aj^ the cofactor of ajk and a, b, c and d are A^-dimensional vectors, a = (oi, 02, oat), 
for example. The following formulas are used frequently in the present analysis [16]: 



an 



yi ... Vn z 



N 



\A\ 



On - i 

am - 1 
-1 



\A\z ^ ] Aji^XjUk, 
j,k=i 



Oijv - 1 1 
-1 1 



{A.9) 



I A(a + b; c + d) 1 = l^(a; c) | + | A(a; d) | + \A{h; c) | + \A{h; d) | , (^.10) 
|A(a,b;c,d)||A| = |A(a;c)||A(b;d)| - |A(a; d)||A(b; c)|, (All) 

N N 

Yl (fj + 9k)ajkAjk = + gj)\A\. {A.12) 

j,k=l j=l 

Formula (A. 11) is Jacobi's identity and formula (A.12) follows from the expansion formulas 
for determinants, XlfeLi ^ikAjk = Sij\A\,Y^k=i '^ki^kj = ^ijl^l- 

Let us now proceed to the proof. First, we modify the determinant \B\. We extract 
a factor 2pj from the jth row of and then extract a factor (1 + Pj)~^ from the jth 
column (j = 1,2, ...,N). Subsequently, the determinant is modified by formula (A.9). We 
extract a factor 1 — pj from the jth row of the resultant determinant and then multiply 
the jth row by a factor 2pj (j = 1, 2, A^). We then find 



g = f,{\A\ + 2|A(-1; q - 1)|) = Ml^l + 2|A(1; 1)| - 2|A(1; q)|). 



(A13) 



where 

q = 



1 '1 ' "■' 1 ) ' 

l-pi I-P2 ^-pnJ 



n7=i(i-p.) 
nf=i(i+p.) 



i = (i,i,...,i), = 



. (.1.14) 



The last hue of (A. 13) is a consequence of formula (A. 10). If we use (A.9), we can rewrite 
(A.13) as 

y = //(|i| + |i(l;l)|-2|i(l;q)|), (A15) 
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where A is a, skew-Hermitian matrix defined by 

^ = {ajk)i<j,k<N, %fe = ie^'Sjk + (A. 16) 

Pj ' Pk 

Since 1^41* = (— 1)^|74|, the comlex conjugate of g becomes 

g* = {-lfl^{\A\ - \A{1; 1)1 + 2|A(q; 1)|). (A17) 

In view of the formulas \A\ = \A\ — \ A{1; 1)| and |A(q; 1)| = \A{q; 1)| which follow from 
(A.9), (A. 17) reduces to 

g* = {-lf^{\A\ + 2\A{q;l)\). (A18) 

Similarly, one has 

r = (-l)^(|A| + 2|A(l;l)|). (A19) 
It follows from (A. 13), (A. 18) and (A. 19) that 

lifg* - fg) = {-1)%[{\A{1- q)| + |A(q; 1)\)\A\ 

+2(Kl;q)| - |A(1;1)|)K1;1)| -2|A||A(1;1)|]. (A20) 

The next step is to calculate the right-hand side of (2.24a). First, applying the differ- 
ential rule of determinant to /, one has 

^ da ^ I ^ f 2 \ 

j,k=l ^ j,k=l j,k=l ^ P^ ~^Pk/ 

By virtue of (A. 8) and (A. 12), we can recast (A. 21) to 

N 

fy^J2P^\A\ + \A{l;p)\, (A22) 
where p — (pi,p2, ■■■,Pn)- A similar calculation leads to 

AT 

9*y = J2p^\B\* + {-If ^l{\A{l■, 1)1 + 1^(1; p)| - \A{q; 1)| - |A(q; p)| - 2|>l(q, 1; p, 1)|). 

{A.23) 
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It follows from (A.18), (A.22) and (A. 23) that 

Dy9*-f = {-If fi[{\A{l; l)|-|A(q; l)|-|A(q; p)|-2|A(q, 1; p, 1)\)\A\-2\A{1; p)||A(q; 1)|] 

{A.2A) 

Using Jacobi's idenity (A. 11) with a = q, b = 1, c = p, d = 1, (A. 24) simplifies to 

Dy9* ■ f = (-1)^/.[(|A(1; 1)1 - \A{q; 1)| - |A(q; p)\)\A\ - 2|A(q; p)||^(l; 1)|]. (A25) 
Referring to (A. 20) and (A. 25), we obtain 

Dy9*-f-l{f9*-rg) 

= {-Ift,{-\A{1; 1)1 + 1^(1; q))| - |A(q;p)|)(|A| + 2|^(1; 1)|). (A26) 

Let P = -\A{1;1)\ + \A{l;q))\ - \A{q;p)\. Applying (A.8) with qj = 1/(1 -p,), P 
becomes 

TV ^ ( \ 

j,k=l j,k=l ^ ^^'^ 

If we extract a factor pi from the /th row of Ajk {I — 1,2, .., N;l ^ j), P is modified as 

(N \ N 

where Aj^ is the cofactor of the {j, k) element of the matrix A defined by 

ie^j 2 

A — {ajk)i<j,k<N, %fe = Sjk H ; • (A. 29) 

Pj Pj+Pk 

Since A is a symmetric matrix, Ajk — A^j. Taking this relation into (A.28), we conclude 
that P — 0. This completes the proof of (2.24a). 

The proof of (2.24b) follows immediately from that of (2.24a) by a symmetry of the 
tau-functions. Indeed, if we exchange the variables y and r and subsequently replace the 
parameter pj by pj^ (j = 1, 2, N), we then see that / from (A. 3) is unchanged whereas 
9 from (A. 4) is transformed to 9*. Thus, under this manipulation, the bilinear equation 
(2.24b) turns out to the bilinear equation (2.24a), which completes the proof of (2.24b). 



29 



Appendix B. An alternative derivation of the 1-soliton solutions 



The 1-soliton solutions take the form of traveling wave 

u = u{X), X = x + cit + xo. (B.l) 

Substituting this expression into equation (1.1) with 1/ — —1 and integrating the resul- 
tant ordinary differential equation once with respect to X under the boundary condition 

u{—oo) = (mod27r), we obtain 

^ (cos-u + c^^ ^ ^ 

Since Ux > 0, we must require that the right-hand side of (B.2) is nonnegative which 
imposes the condition on possible values of Ci. One can see that this condition becomes 
ci > — cos^(ii/2). To proceed, we define a new variable ^ by 

X = y {cosu + ci)d^. {B.3) 

Then, equation (B.2) reduces to 

= ±V(ci + l)2- (cos^/ + Ci)2. (BA) 

In accordance with values of ci, the solutions are classified to several types, which we 
shall now detail. 
1. ci> 

In the case of Ci > 0, (B.4) is integrated through the change of the variable by s = 
tan(ii/2). After an elementary calculation, we obtain 



V ci sinhv^rne 

and 

cosu^- - = 1 5 — ; . (-D-6) 

l + s^ ci sinh^ Vci + U + ci + 1 

Substituting (B.6) into (B.3) and performing the integration with respect to ^, we find 



X ^(^ci + 1)^ - In V^i + 1 + ^^'^^ + U ^ (^ r^^j 

VcTTT - tanh v^cTTT^ 

30 



where is an integration constant. It follows from (B.5) and the boundary condition for 
u that 

u 



= 2 tan~^ ^y^J^^sinh Vci + ^ + ^- (^-^ 



If we put 



ci = i - 1 (0 < pi < 1), ^ = Pii^i - di), rfi = tanhVi, 
Pi 

io = ln('—^)+^ + y,, iB.9) 
\1+Pij Pi 

we can see that (B.7) and (B.8) coincide with (3.4) and (3.3a), respectively. 

2. - 1 < ci < 

In this calculation similar to case 1 gives 



and 



s^±J'-^ } (RIO) 

V — Ci coshyci + l^ 

tan-^ l^^-^cosh v^^TTT^^ + tt, (S.ll) 

X ^{ci + 1)^ - In ^ + + ^ ^^'^^ + ^ 
1 — ^/ci + 1 tanh ^/cl + 1^ 



If we put 



ci = i - 1 (pi > 1), ^ = pi(6 - di), di = tanh ^ — , 

Pi 

Vpi + iy Pi 

we can reproduce the parametric solution (3.3a) and (3.4). 
3. ci ^ 

For the special value ci = 0, integration of (B.3) and (B.4) can be performed readily, 
giving rise to the solution 

u^-+ tan-^(sinh 0, (^-14) 

X = -ln(cosh + Co- (5-15) 
It is easy to confirm that (B.14) and (B.15) coinside with (3.12a) and (3.12b), respectively. 
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Figure captions 



Figure 1. The profile of a regular kink u (solid line) and corresponding profile of v = ux 
(broken line). The parameter pi is set to 0.4 and the parameter yo is chosen such 
that the center position of ux is &t X — 0. 

Figure 2. The profile of a singular kink with the parameter pi — 0.9. 

Figure 3. The profile of a loop soliton with the parameter pi = 2.0. 

Figure 4. The profile of a stationary solution with the parameter pi — 1.0. 

Figure 5. The profile of a kink- kink solution u (solid line) and corresponding profile of v = 

(broken line) for three different times, a: t — 0, h: t — 2, c: t — 3. The parameters 
are chosen as pi — 0.3, p2 — 0.6, = —5, ^20 — 0. 

Figure 6. The profile of a kink-loop soliton solution u for three different times, a: t = 0, b: 
t = 4, c: t = 8. The parameters are chosen as pi = 0.5, p2 = 2.0, ^10 — 0, ^20 — 25. 

Figure 7. The profile of a breather solution v = Ux for three different times, a: i = 0, h: t = 

c: t = 10. The parameters are chosen as pi = 0.2 -f- 0.5 i, P2 = Pi = 0.2 — 0.5 i, ^10 = 

Figure 8. The profile of a sohton-breather solution v = Ux for three different times, a: t — 0, 
h: t — 15, c: t — 25. The parameters are chosen as pi — 0.2 -|- 0.5 i, P2 — P* — 
0.2 - 0.5 i, ps = 0.3, 60 = 60 = 0, 60 = -30. 

Figure 9. The profile of a breather-breather solution v = Ux for three different times, a: t — 0, 
h: t — 40, c: t — 70. The parameters are chosen as pi = 0.1 -|- 0.5i,p2 = 0.16 -|- 
0.8i,p3 =Pt = 0.1 - 0.5 i, P4^P*2^ 0.16 - 0.8i, 60 = Cso = -10, 60 = CIo = 0. 
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Figure 4 
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